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Abstract 

We discuss the generic definition of the r function for the arbitrary 
Hamiltonian system. The different approaches concerning the defor- 
mations of the curves and surfaces are compared. It is shown that 
the Baker-Akhiezer function for the secondary integrable system of 
the Toda lattice type can be identified with the coherent wave func- 
tion of the initial dynamical system. The r function appears to be 
related to the filling of the interior of the classical trajectory by coher- 
ent states. Transition from dispersionless to dispersionful Toda lattice 
corresponds to the quantization of the initial dynamical system. 



1. Recently it was recognized that it is possible to attribute the so called 
r function to any analytic curve in two dimensions [p], p[. This can be done 
as follows. In the simplest case one considers the evolution of the curve when 
the area inside grows linearly while the moments of the curve being fixed. In 
more general situation all moments become new variables ("times") too. It 
appears that in all cases the relevant object which encodes the information on 
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the properties of the deformation of the curve is the quasiclassical r function 
of the Toda lattice with the additional constraint which selects the solution 
to the hierarchy. 

On the other hand it was shown in [[| that it is possible to attribute 
the prepotential which is the quasiclassical version of the r function to any 
holomorphic dynamical system. It plays the role of the generating function 
for the analogue of the S-duality transformation in the context of Hamiltonian 
dynamics. Moreover it was argued that generically there exists a pair of 
dynamical systems such that their coordinates and actions get interchanged. 
Actually one can say that along this viewpoint the quasiclassical r function 
can be attributed to the Riemann surface which is the solution to the equation 
of motion of the holomorphic system. 

The third general setting involving prepotential is the theory of the ef- 
fective actions. It appears that generically effective actions manifest a lot 
of universality property which is governed by the hidden integrability. The 
examples are the identification of the effective action of the N=2 SYM the- 
ory with r function of the quasiclassical hierarchy [[| and similar integrable 
structures behind d=2 theories which have been uncovered in |5], |], [7[. Along 
this way times in the hierarchy are identified with the coupling constants in 
the corresponding theories. In more general terms integrability amounts 
from the deformations of the topological theories or topological sectors in 
nontopological theories by the proper observables 0. 

The goal of this note is to attempt to glue these ideas altogether and for- 
mulate the place of r function for the arbitrary Hamiltonian system without 
any appealing to the integrability. 

2. Let us remind the key points from 0. One can consider complex 
coordinates z, z and the curve determined by the equation 

z = S(z) (1) 

The Schwarz function S(z) is analytic in a domain including the curve. One 
more ingredient to be defined is the map of the exterior of the curve to the 
exterior of the unit disk 

ui(z) = Z - + Y,P 3 z-> (2) 

j 

where uo is defined on the unit circle. It is useful to introduce the moments 
of the curve 

v n =^-jz"S{z)dz (4) 

v = <j) log\z\dz (5) 
which yield the following expansion for the Schwarz function 

S(z) = J2 khz*' 1 + toz- 1 + J2 kv k z- k - 1 (6) 
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Let us define the generating function 

s{z) = d z n{ z ) 



(7) 



with the following expansion 

to(z) = ]T tkz k + hlogz - £ v k z- k k~ l - l/2v (8) 

k=l k=l 

One can easily prove the following relations 

d to n(z) = \ogu(z) (9) 

d tn Sl{z) = (z n (u)) + + l/2(f(u)) (10) 

d t - n Q(z) = + 1/2(S») (11) 

The symbol (S(u>)) + means the truncated Laurent series with only positive 
powers of u kept and the (S(uj)) is the constant term in the series. The 
differential dfl 

dQ = Sdz + logudto + ^^{Hkdtk — Hkdtk) (12) 

provides the Hamiltonians and VL itself can be immediately identified with 
the generating function for the canonical transformation from the pair (z, z) 
to the canonical pair (to,logu). 
The dynamical equations read 

d tn S(z) = d z H n (z) (13) 

d fn S(z) = d z H n (z) (14) 

and the consistency of (13) (14) yields the zero-curvature condition which 
amounts to the equations of the dispersionless Toda lattice hierarchy. The 
first equation of the hierarchy reads as follows 

dl fl <f> = d to e d *ot> (15) 

where d to (p = 2logr. The Lax operator L coincides with z{uj) and its eigen- 
function - Baker- Akhiezer function looks as follows ^ = eft. Hamiltonians 
are expressed in terms of the Lax operator 

H k = {L k ) + + l/2(L fe ) (16) 

3. In this section we shall interpret the evolution above in terms of 
dynamics on the phase space assuming that z, z pair yields the phase space 
of some dynamical system. The curve itself corresponds to the energy level 
of this dynamical system. With this setup it is clear that Poisson bracket 
between z and z is fixed by the standard symplectic form . The next step 
is the identification of logw as angle variable hence the identification of the 
area inside the curve t as the action variable is evident. 
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The holomorphic variables on the phase space are usually considered as 
the operators in the Fock space so we identify z = b;z = b + . Typically 
the creation- annihilation operators are expressed it terms of the (p, q) vari- 
ables however for our needs we assume the action-angle representation. The 
variable z = b can be expressed in terms of action-angle variables in the par- 
ticular dynamical system and just this expression provides the Toda lattice 
Lax operator 

z(logw, t ) = L(logw, t ) (17) 
Baker- Akhiezer (BA) function is solution to the equation 

Ltf(z,to) = «tf (18) 

and the complex z plane is identified with the surface where the spectral 
parameter of the Toda system lives on. 

Now we turn to clarification of the meaning of the BA function in the 
generic dynamical system. The answer follows from the equation above; BA 
function is nothing but the coherent wave function in the action representa- 
tion. Indeed the coherent wave function is the eigenfunction of the creation 
operator 

m = b^ (19) 

The main properties of the coherent states can be summarized as follows 
The coherent state can be represented in the form 

\a>=D(a)\0> (20) 

where |0 > is the normalized vacuum vector annihilated by the operator b 
and the D(a) = exp(ab + — ab) yields the action of the Heisenberg-Weyl 
group. The operators D(a) obey the following relation 

D(a)D((3) = e 2iIm{a ^D((3)D(a) (21) 

The coherent states \a > and \/3 > are not orthogonal moreover the system of 
the coherent states [\a >] is overcompleted generically. To get the complete 
system one can proceed as follows. The lattice on the complex plane a nm = 
nuii + vniOi can be introduced with the area of the elementary cell equals to 
the elementary Planck cell ir. In this case the system of the coherent states 
\ot nm > can be expressed in terms of the standard theta-functions on the 
torus and turns out to be complete. 

Note that the system of the coherent states can be defined for the generic 
group as follows 

\<iP g >=T(g)\ip > (22) 

If H is the stabilizer of the vector \ip > the coherent state is defined by the 
point x=x(g) from the homogeneous space G/H corresponding to the element 
g. Moreover the state corresponding to \x > can be identified with the one 
dimensional projector P x = \x >< x\. 

Now we are ready to recognize the meaning of the generating function 
Q. From the equations above it becomes clear that it is the generating 
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function for the canonical transformations from the b, b + representation to 
the angle-action variables. Let us also remark that the role of the Orlov- 
Shulman operator becomes transparent. This operator looks like b + b and its 
eigenvalue counts the "number of particles". 

Having identified the BA function for the generic system let us consider 
the role of the r function in the generic case. To this aim it is convenient to 
use the following formulae for the r function 

T (t, W) =< t,t\W > (23) 

where the bra vector depends on times while the ket vector is fixed by the 
so called point of Grassmanian. One more useful representation is provided 
by the fermionic language 

T(tiW) = <mn^M\w> (24) 

where A(z) is Vandermonde determinant. 

Let us briefly comment on the definition of the point of the Grassmanian 
W. Generally speaking the Grassmanian itself can be considered as a collec- 
tion of all fermionic Bogolyubov transforms of the vacuum state. Hence we 
can say that W belongs to the Grassmanian if it is annihilated by some linear 
combination of the fermionic creation and annihilation operators. Equiva- 
lently one could consider the following definition 

\W >= S\0 >,S = exp^Amd-n-xpii-m-ifr (25) 

nm 

that is S can be considered as the element of GL(oo, C). 

The consideration above suggests the following picture behind the defi- 
nition of the r function. The peculiar classical trajectory of the dynamical 
system yields the curve on the phase space. Then the domain inside the tra- 
jectory is filled by the coherent states for this particular system. Since the 
coherent state occupies the minimal cell of the phase space the number of the 
coherent states packed inside the domain is finite and equals N. Since there 
is only one coherent state per cell for the complete set it actually behaves 
like a fermion implying a kind of the fermionic representation. 

Therefore we can develop the second dynamical system of the Toda lattice 
type based on the generic dynamical system. The number of the independent 
time variables in the Toda system amounts from the independent parameters 
in the potential in the initial system plus additional time attributed to the 
action variable. The situation is similar to the consideration of the pair of 
the dual dynamical systems in ||. The subtle point in that paper is the 
choice of the Hamiltonians for the dual systems. It seems that consideration 
in this note provides the unified viewpoint on this issue. Indeed formulae 
above suggest the set of commuting Hamiltonians for the dual system. Let us 
emphasize that the choice of the particular initial dynamical system amounts 
to the choice of the particular solution to the Toda lattice hierarchy. 
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It is important that the energy level of the dynamical system doesn't 
develop the discontinuity in some range of the energies. This happens only if 
the system approaches separatrices generically existing on the phase space. 

Consider the example of the oscillator. Corresponding phase trajectories 

E = p 2 + q 2 (26) 

are circles on the phase space (ellipses in the generic case) and the natural 
complex variables are 

z = p + iq (27) 

If we calculate the expression for vq the following prepotential for the oscil- 
lator emerges 

T = l/2I 2 logI -3/4I 2 (28) 

where I is the action variable. To get the first dispersionless Toda equation 
one has to consider the phase trajectories with the shifted center. The posi- 
tion of the center is parameterized by times £i, t\. If the Toda equation with 
three nonvanishing times is solved the prepotential above can be reproduced. 
Let us also note that the dual Hamiltonian in the sense of [[3] is 

H = b = Vhe l9t ° (29) 

The wave function of the dual system coincides with the coherent wave func- 
tion in the action representation. 

The prepotential for the complex system is defined by its spectral curve 

dT , s 

a D = (30) 

where and a are the integrals of the meromorphic differential over A 
and B cycles. The variable a is just the action variable. The S duality 
transformation a — > actually maps the region of small energies to the 
region of the large energies. Prepotential obeys the Matone type || relation 
which has the meaning of the Ward identity 

dT 

constH (31) 



logA 

where H is the Hamiltonian of the complex system and A is the scale factor. 
The simple calculation shows that such relation doesn't 't hold in the real 
case. 

Let us remind how the prepotential provides the S-duality transformation 
for the complex dynamical system |J. The action variables in dynamical 
system are the integrals of meromorphic differential A over the A-cycles on 
the spectral curve. The reason for the 5-cycles to be discarded is simply the 
fact that the S-periods of A are not independent of the A-periods. On the 
other hand, one can choose as the independent periods the integrals of A over 
any lagrangian subspace in Hi(Tb] Z). 

This leads to the following structure of the action variables in the holo- 
morphic setting. Locally over a disc in B one chooses a basis in H\ of the 
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fiber together with the set of A-cycles. This choice may differ over another 
disc. Over the intersection of these discs one has a Sp(2m, Z) transformation 
relating the bases which has a natural interpretation as a S-duality transfor- 
mation. Altogether they form an Sp(2m, Z) bundle. It can be easily shown 
the two form: 

dP A dlP (32) 
vanishes. Therefore one can always locally find a function T , such that: 

IF = % (33) 

The angle variables are uniquely reconstructed once the action variables are 
known. Since the real system to some extend can be considered as a real 
section of the dynamical system one could expect the analogue of the S- 
duality for the real system too. 

An important outcome from our identification is the natural recipe for 
the transition from the dispersionless to dispersionful Toda lattices. Indeed 
at the first step of quantization the system one imposes the Bohr-Sommerfeld 
quantization condition on the action variable. This is in agreement with the 
discretization of to variable in Toda lattice. The Planck constant in the initial 
system plays the transparent role in the Toda lattice. 

4. Let us argue now that there is clear correspondence between the inte- 
grable dynamics behind the deformation of the curves which are the phase 
space trajectories for the real dynamical systems and the deformation of the 
surfaces which represent the trajectories of the complex dynamical systems. 
Remind that to formulate the integrable dynamics behind the deformation 
of the curves one introduces the following main objects; generating function 
Q, t function depending on time variables t k which are the moments of the 
exterior of the classical trajectory as well as the dual variables v k representing 
the moments of interior. 

It is not a difficult task to recognize their counterparts in the case of the 
surfaces. First let us represent the time variables in the following form 

v k = £z k dn(z) (34) 

t k = I z- k dVt(z) (35) 
and compare it with the definition of times in the case of surfaces 

T„ = res (= „CdSsw (36) 

^ = 2L res °*-" ds ™ (37) 

where Ssw in the context of the N=2 SYM theories is the so-called Seiberg- 
Witten meromorphic differential whose derivatives with respect of the moduli 
of the curve are the holomorphic differentials. We see that the variables in 
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cases of curves and surfaces are differed just by the Legendre transform. The 
variable to is nothing but the variable 
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dS (38) 



where the integration is over the A cycle on the spectral surface is the com- 
plex dynamical system with only one degree of freedom is considered. The 
definition of conjugated variables v n is more delicate issue since naively there 
are no the direct analogue of the integrals over the B cycles 



a D = <j> dS (39) 

in the real case. However it can be seen that the analogous variable for the 
one degree of freedom can be identified with variable vq from Jl|, . If there 
are several degrees of freedom the corresponding variables are defined as 
the integrals over the complete set of A cycles. 

It is interesting to recognize the real analogue of the modulus of the 
spectral curve in the complex case 

r = (40) 
da 

which plays the role of the effective coupling constant in the field theory 
context. In the problem under consideration the corresponding variable is 

logr = P (41) 
y dt y 1 

and also exhibits the behaviour similar to the perturbative regime in the field 
theory. In what follows we shall see more indications that dynamics of the 
curve has a pertubative regime as a field theory counterpart. 

In the context of the Seiberg-Witten solution to the N=2 SUSY YM the- 
ory the Whitham dynamics has a transparent renormalization group meaning 
||10|| . Indeed the derivative of the prepotential with respect to the first time 
in SUSY YM context looks as follows 

^ = a^- -IT = 2H(a) (42) 
0Ti da K 1 y ' 

where the term at the r.h.s. coincides with the Hamiltonian for the dy- 
namical system providing the corresponding spectral curve. The dependence 
H(a) on the action variable is the intrinsic characteristics of the dynamical 
system, moreover this term has an anomalous nature from the field theory 
viewpoint. The analogous equation has been also found for the real case [J] 
where the observed quadratic dependence on the action in the anomaly term 
corresponds to the perturbative limit in the field theory framework. 

Given such RG interpretation in the complex case one could look for the 
reasonable RG structure behind the real dynamical system. The following 
starting point can be chosen. Assume that the RG scale is fixed by the energy 
corresponding to the classical trajectory or by the variable to = A*. Then the 
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phase space is naturally divided into the "IR theory" inside the trajectory 
and the "UV theory" outside. Therefore a kind of the RG problem can 
be formulated. To this aim introduce the classical expectation value of the 
observable in the theory 

<0>= J dpdqO(p, q)p(p, q) (43) 

where p(p,q) is the phase space density (Wigner function). The vacuum 
expectation value at the scale p is naturally introduced if the integration 
over the phase space up to the energy p is performed. 

Now we can formulate the following observation - correlators in UV theory 
in the outer region amount to the coupling constants in the inner region and 
vice versa 

r(t , t k ) =< exp t k O k > (44) 

where the matrix element is taken in the IR theory. The key point is that 
usually effective action S(p) identified with the prepotential is defined via 
integrating out all modes up to the scale p. From U it becomes clear that 
is has to be done in a clever way; one has to integrate out the part of the 
phase space outside of the phase trajectory corresponding to the energy E. 
Then the Riemann-Hilbert problem is defined on the RG scale p - one has to 
divide the whole theory into UV and IR parts consistently with the RG flows. 
Remark that interpretation of the area inside the curve as the RG parameter 
agrees with the interpretation of the size of the matrix in the matrix model 



as RG scale |12|]. Let us emphasize that the picture emerged is formally 
identical to the consideration in jTD| where r function approach amounts to 
the explicit calculation of the RG dependence of the correlators. Along this 
viewpoint one could expect that the similar correlators which now carry the 
information about the spectrum of the system could be found. Therefore 
instead of the determination of the RG invariant characteristics of the field 
theories we expect the RG like evaluation of the spectral invariants which 
would distinguish different behaviour of the Hamiltonian systems. 

Let us emphasize that the Toda lattice equations are written for the 
variable b(uj) which yields the dependence of the creation operator on the 
angle variable. If only zero time is involved then the evolution doesn't change 
the system however with all times switched on the dynamics proceeds in the 
space of the dynamical systems or can be interpreted as the RG flows. The 
integrals of motion in the Toda system are the "classical correlators" in the 
IR theory. 

5. It was shown in |l| that the r function of the curve allows two ad- 
ditional complementary interpretations. First, it has the interpretation of 
the energy of the fermions with the unit charge inside the curve with the 
Coulomb interaction taken into account. Secondly it can be considered as 
the partition function of the normal matrix model with the potential yielding 
the set of times in the Toda lattice. 

Let us speculate on the appearance of the matrix model in a given context. 
To this aim consider the phase space as a two dimensional noncommutative 
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plane. Evidently the Planck constant h fixes the scale of noncommutativity 
due to the canonical commutation relation 



[xi, x 2 ] = ih (45) 

The Moyal multiplication law for the functions on the phase space reads 

f(x) * g(x) = exp(ihd x d y ) x=y g(y)f(x) (46) 

Now let us remind the notion of the Morita equivalence. Qualitatively it 
means that for the rational noncommutativity (h = p/q) the noncommu- 
tative theory is equivalent to the commutative theory on the rescaled non- 
commutative manifold. The functions on the commutative manifold become 
matrix valued with the size of the matrix fixed by the noncommutativity. 

Let us apply these circle of ideas to the generic dynamical system. Sup- 
pose that the Planck constant is a rational number h = ■h. If N is large 
we are dealing with the quasiclassical approximation. Applying the Morita 
equivalence, all functions become N x N matrix valued and the system be- 
comes effectively classical. Therefore Morita transform maps the quantum 
system into the classical matrix theory. The " classical" partition function of 
the matrix model is the prepotential for the dynamical system. 

One more indication of the effective classicality of the matrix model fol- 
lows from the matrix relation [M + ,M] = which corresponds to the "clas- 
sical" commutation relation. On the other hand the corresponding Planck 
constant in the Morita equivalent system is h = N which indicates the regime 
of the strong noncommutativity on the phase plane. Along this viewpoint 
one could expect the appearance of the peculiar objects relevant for the such 
regime, namely the noncommutative solitons |13[ which in the simplest case 



correspond to the projector operators \n >< n\. One the other hand we 
have already remarked that the coherent states also have the description in 
terms of the projector operators on the homogeneous spaces. This implies 
the intriguing conjecture to consider the matrix model as the system of the 
interacting noncommutative solitons. 

Let us make one more remark concerning the different interpretations of 
the prepotential. It has a interpretation of the energy in the matrix model, 
action for the c=l string and somewhat involved interpretation in the context 
of particle dynamics. This could presumable have the following picture be- 
hind; it is energy for the membrane and therefore the action for string. With 
this picture in mind the matrix model corresponds to the M(atrix) model 
representation for the membrane. This is in a rough agreement with the 
chain of T duality transformations of the phase space. Let us also note that 
in the context of c=l string "times" correspond to the vacuum expectation 
values of tachyons with the different momenta which could be related to the 
known phenomena of an appearance of the noncommutative solitons(coherent 



states) as the result of the tachyon condensation [I| . 

6. To conclude we have argued that the Toda lattice governs both the de- 
formations of the complex surfaces and analytic curves. These deformations 
can be considered as the motion in the space of the complex and real dynam- 
ical systems respectively or as a kind of RG flows in a peculiar system. The 
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dynamics with respect to the time to can be considered as the RG evolution 
with to playing the role of the scale factor. One could also expect that the 
discrete symmetries could provide an important tool for investigation of the 
RG flows like in [|TjJ. 

The Baker- Akhiezer function was definitely identified as the coherent 
wave function in the action-angle representation and the generating func- 
tional has the interpretation of the generating function of the canonical trans- 
formation from holomorphic to action-angle representation. The role of the 
quantization of the initial dynamical system in the integrability framework 
is clarified. 
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fellowship of CNRS 2000. 
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